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Abstract: We study two possible explanations for short baseline neutrino oscillation anoma- 
lies, such as the LSND and MiniBooNE anti-neutrino data, and for the reactor anomaly. The 
first scenario is the mini-seesaw mechanism with two eV-scale sterile neutrinos. We present 
both analytic formulas and numerical results showing that this scenario could account for the 
short baseline and reactor anomalies and is consistent with the observed masses and mixings 
of the three active neutrinos. We also show that this scenario could arise naturally from 
an effective theory containing a TeV-scale VEV, which could be related to other TeV-scale 
physics. The minimal version of the mini-seesaw relates the active-sterile mixings to five real 
parameters and favors an inverted hierarchy. It has the interesting property that the effec- 
tive Majorana mass for neutrinoless double beta decay vanishes, while the effective masses 
relevant to tritium beta decay and to cosmology are respectively around 0.2 and 2.4 eV. The 
second scenario contains only one eV-scale sterile neutrino but with an effective non-unitary 
mixing matrix between the light sterile and active neutrinos. We find that though this may 
explain the anomalies, if the non-unitarity originates from a heavy sterile neutrino with a large 
(fine-tuned) mixing angle, this scenario is highly constrained by cosmological and laboratory 
observations. 
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1. Review of the Neutrino Oscillation Anomalies 

Neutrino oscillation experiments with solar, atmospheric, reactor, and accelerator neutrinos 
have established a standard picture of three-flavor neutrino mixing and masses [1] : 

7.05 x l(T 5 eV 2 < Am^ < 8.34 x 1(T 5 eV 2 , 
2.07 x 10" 3 eV 2 < lAmjy < 2.75 x 10~ 3 eV 2 , 
0.25 < sin 2 0i2 < 0.37, 0.36 < sin 2 <9 23 < 0.67, sin 2 13 < 0.056 at 90% C.L. 

However, over the past ten years, several anomalies have been observed in the short baseline 
experiments, suggesting deviations from the standard three-flavor picture. They are 

• LSND [2]: the search for 9u — > 9 e conducted using the Liquid Scintillator Neutrino 
Detector with L/£'(m/MeV) = 0.4 - 1.5 (m/MeV), where L is the distance between 
the neutrino source and detector and E is the neutrino energy. A 3.8(7 excess of 9 e 
candidate events was reported. The 9^ arises from 7r + /^ + decay at rest: tt + — > [J^Vp, 
and [i^ — > e + v e 9 fjL . The signal is defined as a positron with a 52 MeV endpoint from 
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v e p — > e + n with a correlated 2.2 MeV photon from neutron capture on a free proton 
np —> drf. Another experiment, KARMEN [3], also searched for &„ — > v e using from 
[i + decaying at rest, with L/E(m/MeV) ~ 0.4 - 1.8 (m/MeV). The result is consistent 
with zero oscillation events and leads to a 90% C.L. upper limit on the oscillation 
probability Pp -+p e < 6.5 x 10~ 4 . There is a small region of parameters consistent with 
both LSND and KARMEN [4]. 

• MiniBooNE [5, 6, 7]: searches for both z/„ — > u e and — > v e were carried out with 
L/E = 0.2 - 2.6 (m/MeV). In the neutrino search [5], no excess of v e was found above 
475 MeV 1 . In the anti-neutrino search reported in 2010 [6], an excess of v e events above 
475 MeV was observed with 2.75<r. The 2011 updated result represents an increase in 
statistics of 52% and still favors LSND-like — > v e oscillations over a background only 
hypothesis at 91.1% C.L. in the range L/E = 0.2 - 1.13 (m/MeV) [7]. 

• Reactor anomaly [8]: A re-evaluation of the expected mean reactor anti- neutrino flux 
yields a prediction 3% higher than what was previously assumed. This implies that 
all reactor neutrino experiments with L < 100m have observed a deficit of T> e events 
compared to the theory prediction, of 5.7% on average, at 98.6% C.L. 

These three experiments could potentially be interpreted as oscillations between eV-scale 
sterile neutrinos v s and the active neutrinos 2 . While the 3+1 scheme (with only one eV-scale 
u s ) could still not be consistent with all data, Kopp, Maltoni and Schwetz (KMS) [10] and in- 
dependently Giunti and Laveder (GL) [11], showed that with the new reactor flux evaluation, 
the global fit of these and other short baseline experiments improves greatly in a 3+2 scheme 
with two eV-scale sterile neutrinos. This is because the three active neutrinos are essentially 
degenerate on the relevant eV-scale, requiring two sterile neutrinos for a CP-violating dif- 
ference between za, — > v e and — > v e oscillations. As already pointed out in [12], the 3+2 
scheme could be accommodated in a mini-seesaw paradigm, which we will prove in more detail 
in this paper. Moreover, we will show that the mini-seesaw paradigm could easily be realized 
in an effective theory involving a scalar with a TeV-scale vacuum expectation value (VEV). 
The most serious constraints on this scheme are from cosmological observations. Current 
data favors one or two light sterile neutrinos, but their masses are strongly constrained: one 
eV-scale sterile neutrino may be barely allowed, but two are inconsistent with all observations 
unless certain cosmological assumptions are relaxed [13]. 

Another possible explanation for the short baseline anomalies is a model in between the 
3+1 and 3+2 schemes, in which one sterile neutrino is light and the other is heavy [14]. In this 
case, only one oscillation length is comparable to the range of the experiments. The heavy 
neutrino is either not kinematically produced or its much shorter oscillation length is averaged 
over. The situation could be described by a non-unitary mixing matrix [15, 16] for the active 

1 A 3cr excess of electron-like events at lower energy, apparently unrelated to the LSND signal, was reported. 
The origin of the excess is unknown. 

2 Alternative explanations are summarized in [9]. 
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neutrinos and one light sterile neutrino. Without the unitarity constraint, CP violation could 
be present even for two- neutrino oscillations [17]. Naively one might expect cosmology to be 
more friendly to such a scheme as only one eV-scale sterile neutrino is present. However, we 
found that the heavy sterile neutrino still needs to have a relatively large (fine-tuned) mixing 
angle with the active ones to generate large enough CP violation if one wants to explain all the 
data. Thus, this scheme is also highly constrained by cosmological/astrophysical observations 
as well as by laboratory experiments such as meson decays. 

Finally we mention another anomaly that receives relatively less attention, the "Gallium 
anomaly". The SAGE [18] and GALLEX [19] experiments used megacurie sources of 51 Cr 
and 37 Ar to calibrate v e + 71 Ga — > 71 Ge + e~ . The data from these two experiments are 
consistent with each other and show a deficit compared to the theory predictions, which can 
be interpreted as a 2.7a indication of v e — > u s , where u s is a sterile neutrino with Am 2 ~ 2.24 
eV 2 [20]. If true, this has a subtle implication for the estimate of v e background in the 
MiniBooNE neutrino oscillation channel. However, as shown by [21], the KARMEN and 
LSND cross section measurements for v e + 12 C — > 12 N gs + e~ exclude the best fit point for the 
Gallium anomaly at 3.6cr. Thus, we do not put this anomaly on our motivation list. 

The paper is organized as follows. In Section 2, we first review the basic formalism for 
neutrino oscillations involving sterile neutrinos and the global best fits from the literature. 
Then we examine the mini-seesaw mechanism and propose an effective theory that realizes 
the mini-seesaw naturally. A minimal version of the mini-seesaw, in which the active neutrino 
masses are due entirely to mixing, is especially predictive and is examined in detail. It strongly 
favors an inverted hierarchy (IH) but implies a vanishing effective neutrinoless double beta 
decay mass mpp. It is shown that this minimal version can accomodate the observed active 
neutrino mass differences and mixings and can successfully fit the LSND and MiniBooNE 
data. We conclude the section by discussing the cosmological constraints and possible ways 
to evade them. In Section 3, we discuss another possible scenario that may explain the short 
baseline anomalies with one light and one heavy sterile neutrino. We comment on the rather 
stringent cosmological and astrophysical constraints on this scenario and possible evasions. 
We conclude in Section 4. Further details of the minimal mini-seesaw fit are given in the 
Appendix. 

2. The Mini- Seesaw Mechanism 
2.1 Global Fits 

We first review the basic formalism for neutrino oscillations. Neglecting the masses in the ac- 
tive neutrino sector and assuming that only two sterile neutrinos participate in the oscillation, 
the probability of electron neutrino appearance in a muon neutrino beam is 

Pv^v e = 4|[/ e4 | 2 |£/ M 4| 2 sin 2 X4 + 4|C/ e5 | 2 |C/ / ,5| 2 sin 2 X5 

+8|J7 e4 ||?7 / ,4||[/ e 5||£/ M 5| sina;4sinx5Cos(x5 - x 4 - 5), (2.1) 
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where U a i are the mixing matrix elements between active neutrinos v a (a = e, fx) and steriles 
V{ (i = 4, 5), and the phases X{ are 



1.27 



m i 
eV2 



L/E \ 
m/MeV J 



(2.2) 



In propagation through matter the phases may also be altered by forward scattering due to 
the weak interactions and thus are also sensitive to E. The matter effects are generically 
small in short baseline experiments unless exotic forces exist, which we will neglect in this 
paper. The physical CP-violating angle is defined as 



8 = are 



u e4 u, 



(2.3) 



For anti-neutrino oscillations, one only needs to replace 5 by —5 in Eq. 2.1. 

As shown by [10] and [11], with the new reactor flux prediction, the global fit to all 
oscillation data (based on 2010 MiniBooNE anti-neutrino data) improves considerably when 
the existence of two light sterile neutrinos is assumed. The best fit points, given in Table 1, 
prefer two steriles with masses of order 1 eV; the active-sterile mixing angles of order 0.1; 
and a non-zero CP violating phase 5 in the sterile sector. 







Am 2 51 


\u e 4 


\u e5 \ 


l^4| 


\u„ 5 \ 


S/tt 


X 2 /dof 


KMS 
GL 


0.47 
0.90 


0.87 
1.60 


0.128 
0.13 


0.138 
0.13 


0.165 
0.13 


0.148 
0.08 


1.64 
1.52 


110.1/130 
22.2/5 



Table 1: Best global fit points of KMS [10] and of GL [11]. (The GL result is obtained from parameter 
goodness of fit where the number of degrees of freedom corresponds to the number of parameters in 
common to the data sets. For more details, see [22].) 



2.2 Active-Sterile Neutrino Mixing 

Most extensions of the standard model (SM) which allow nonzero neutrino mass 3 involve 
sterile (i.e., <S?7(2)-singlet) neutrinos. Consider, for example, the case of three left-chiral active 
neutrinos i/£ and n left-chiral sterile antineutrinos N® c , and their respective CP conjugates 
and Nft, where the superscript indicates weak eigenstate. In general, these will have a 
mass term 

-H(^K)(|)> 

where Mt = M^ and Ms = Mj are 3x3 and n x n-dimensional symmetric Majorana mass 
matrices, while Mrj is a 3 x n-dimensional Dirac mass matrix. Mt, Md, and Ms violate weak 
isospin by 1,1/2, and units, respectively. In the pure Majorana (Md = 0) and pure Dirac 



3 For general reviews, see, e.g., [23, 24, 25, 26]. 
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(Mt = Ms = 0) cases there is no mixing between v® L and N® c (or between z/^ c and N^). In 
the ordinary Type I seesaw model (with the Ms eigenvalues much larger than Mb and Mt, 
typically > O(TeV)) there is mixing, but it is very small. Moreover, the heavy, predominantly 
sterile, states decouple from the low energy theory. Therefore, mixing between active and 
light sterile neutrinos of the same helicity, which could account for the anomalies described in 
Section 1, requires the simultaneous presence of small Dirac mass terms and small Majorana 
mass terms (see, e.g., [27]). We will therefore consider an effective low energy theory in 
which all of the entries in (2.4) are very small (usually < O(eV)) after integrating out any 
heavy states. Especially attractive is the mini-seesaw, which is similar to the ordinary seesaw 
except that the eigenvalues of Ms are in the eV range (and much larger than the elements 
of Md and Mt)- This typically leads to relations between the active-sterile mixings and the 
mass eigenvalues consistent with the observations. For example, in the one- family version of 
(2.4), the 2x2 mass matrix is 

M= (m T m D \ 
ym D m s J 

Taking mj- = for simplicity and ms S> mp, the physical masses of the predominantly 
active and sterile states are m\ ~ m D /ms and m.2 ~ mj, respectively, with an active-sterile 
mixing 6 ~ mrj/ms ~ (m-i/m^) 1 / 2 . Thus, = 0(0.1 eV) and ms = 0(1 eV) implies 
mi = 0(0.01 eV) and 9 = O(0. 1). The extension to several families and the roles of family 
mixing and Mt are discussed below. 

2.3 Mini-Seesaw from Higher-Dimensional Operators 

Barring fine-tuning, the small values of Mp and Mt,s needed for active-sterile mixing or for 
the mini-seesaw most likely imply that they are suppressed by additional gauge, global, or 
discrete symmetries compared to the simplest expectations. The remaining small elements 
may usually be described 5 by higher-dimensional operators (HDO) involving powers of S/A, 
where S represents a SM singlet field (or fields) whose VEV breaks the symmetry. A is 
the new physics scale, typically involving heavy particles or string excitations that have 
been integrated out. The lowest-dimensional operators 6 yielding a realistic spectrum are (in 
superpotential notation) 

(H L) 2 9 S 2 

W T ~ " , W D ~ —H U LN, W s ~ ^N 2 . (2.5) 



4 It could instead be due to two distinct kinds of small Dirac mass terms, such as one which links active 
and sterile neutrinos, and another which links sterile left and right-chiral fields. 

Another possibility is that the mass terms are exponentially suppressed by nonpertubative effects, such as 
D-brane instantons. For reviews, see [28, 29]. 

6 The operator approach for neutrino masses has been developed generally in [27, 29, 30, 31]. Specific models 
that can lead to such operators include string or U(l)' motivated models [27, 32], mirror worlds [33, 34], gauge- 
mediated supersymmetry breaking [35], compositeness [36], dynamical electroweak symmetry breaking [37], 
warped extra dimensions [38], and variations on conventional seesaw and flavor models [39, 40, 4f, 42, 43]. 
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In the mini-seesaw limit this implies 

( M T M D \ ( O(^) O(^) \ ( 0(0.01) 0(0.1) 



\Ml M s ) \o(f)0(^)J \ 0(0.1) 0(1 

where v ~ 246 GeV is the electroweak scale and we have taken (S) = O(TeV) and A ~ 10 15 
GeV. Note that if both Wd and Ws are allowed by a multiplicative symmetry then so is Wt, 
and that in the mini-seesaw limit the mixing-induced contribution to the light eigenvalues, 
Mp/Mg ~ ^ 2 /A, is generically comparable to that from Wt- 

2.4 Active-Sterile Mixing Parameters in the Mini-Seesaw 

The implications of the general mass terms in (2.4) and its (mini-)seesaw limit have been 
discussed in detail in [44, 12, 45, 46, 47, 48]. Here, we recount the relevant features for our 
analysis. In the general case the symmetric mass matrix in (2.4) can be diagonalized to yield 
3 + n Majorana mass eigenstates 

vl 

where A U L is a (3 + n) x (3 + n) unitary matrix and vl is a (3 + n)-component vector. The 
analogous transformation for the R fields is A V R = A 1 ^ . In the important special case My = 
and n < 3, there will be 3 — n massless active neutrinos 7 . Thus, n = 1 or 2 sterile neutrinos 
and Mt = would be candidates for describing the NH or the IH for the active neutrinos, 
respectively. 

In the seesaw limit, where the eigenvalues of M$ are both large compared to all entries 
of and Mp, one has 

A« = A'?=( A ? iW, 






where (A^ ) are 3 x 3 (n x n) unitary matrices and 

B *t = [ 1 -M D M~ S 
L \ M~ 1] M ] D I 



l 



One finds that 



B 



v] / M T M d ) B u* = ( M t - M D M^M T D 
1 Ml M s i L 1 M s 



Without loss of generality, we can choose a basis in which Ms is diagonal, with real and 
positive elements. We will first consider the limit Mt = 0, which we refer to as the minimal 
mini-seesaw. The mass matrix for the three active neutrinos is then 



MdM^M^ = DD T , (2.6) 



7 There can be additional massless states for singular Mo, but we will not consider that case. 
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— 1/2 

where D = MdM s is a (3 x n)-dimensional complex matrix, and the extra minus sign has 
been absorbed in a redefinition of the weak eigenstate phases. m u can also be expressed as 



A v L m d Af = LL 



(2.7) 



where is a diagonal matrix of the three light neutrino mass eigenvalues m», which can 
be taken to be real and positive by appropriate choices of mass eigenstate phases, and L 
A v L mJ . It is convenient to view L as (3 x n)-dimensional: for n < 3 one can simply remove 
the 3 — n columns of zeros corresponding to rrii = 0, while for n > 3 one can add n — 3 extra 
columns of zeros. We work in a basis for which the charged lepton mass matrix is diagonal, 
in which case A V L is just the active neutrino (PMNS) matrix [49, 50] 



1 

C 2 3 S23 
-S 2 3 C 2 3 






( 


C12 


•S12 




) 




-S12 


C12 






{ 











(2.8) 



where Cj,- = cos #™ and Sj,- = sinfe are leptonic mixings, p is the CP-violating phase, and 



$ = diag(e* Ql , e ia2 , e' tQf2 ) is a diagonal phase matrix. The Majorana phases cti,i = 1,2,3 are 
not observable in ordinary 3-flavor oscillations, but they do affect the active-sterile mixing. 
The charged lepton phases can be chosen so that one of the at is zero. We will use this 
freedom to choose a\ = 0. Then, as is unobservable for 7713 = 0. The observational data is 
consistent with the tri-bimaximal form s\2 = l/v^, #13 = 0, #23 = 7r /4 [51, 52], i.e., 



( 



\ 







\ 



2/ 



( 

V 



0.82 0.58 
-0.41 0.58 0.71 
0.41 -0.58 0.71 



<I>, 



(2.9) 



except possibly for recent hints from T2K, MINOS, and Double CHOOZ for a nonzero #13 [53, 
54, 55]. We will use this form in our numerical illustrations, but have verified that small 
modifications are not important for our purposes. 

The active-sterile neutrino mixing is described by the 3 x n matrix 



U 



iMoMe 1 = iDMo 1/2 , 



where the i results from the phase redefinitions mentioned above. From (2.6) and (2.7), the 
most general solution for D in terms of L is 



D = LR(zk), 



(2.10) 



where R is an orthogonal n x n complex matrix, which depends on n(n — l)/2 complex 
parameters as well as signs. Thus, the mixings can be predicted in terms of the PMNS 
matrix, the light and heavy mass eigenvalues, and the z^. 
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For n = 1 sterile neutrino one has a NH, with mi = m2 = 0, and with 7773 ~ 0.047 eV 
determined from the atmospheric and long baseline neutrino oscillations. The active-sterile 
mixings are given by 

U a4 = i— ^ = ±iA v L a3 ,/— ^, a = e,fx,T. (2.11) 



M 4 u V M 4 

Because of the stringent upper limit on A 1 ^ 3 < 0.24, \J e4 is too small to account for the 
LSND/MiniBooNE anomaly [12] (even ignoring the need for CP violation to differentiate 
from and the lack of a Solar mass splitting between mi and m^). 
For n > 1 



Urvi ~~ ^ . , == %A T ^ \/7Tl7 H>ji , 

Mi v/M 



(2.12) 



where a = e, /i, r; j = l,-- ,n; and i = 4, • • • , 3 + n. In particular, for n = 2 there is 
one massless neutrino. In principle one could have either a NH or IH for the three light 
neutrinos. However, for the normal case it is difficult to obtain large enough active-sterile 
mixings because of the small values of Aj^ 3 and mi- We therefore consider an IH with 771,3 = 0, 
with the active-sterile mixings determined by the masses and PMNS matrix, as well as one 
complex parameter 8 z, i.e., 

R(z)=( C ° SZ S[nZ ). (2.13) 
\ — sin z cos z I 

Because of the possible Majorana phases in A V L and because R is complex, there is room 
for CP violation in the active-sterile mixing. Note that all values of the PMNS parameters 
and of the mass eigenvalues (with 7773 = 0) are consistent with the minimal mini-seesaw 
parametrization. 

To test whether Eq. 2.12 could account for the LSND and MiniBooNE results, we take 
the experimental data to construct a x 2 function 

^(^ theory (.,a 2 ,M 4 ,M 5 )-^) 2 
X 2 (z,a 2 ,M A ,M 5 ) = ^2^ ^ J -, (2.14) 



"1 



where P° xp is the oscillation probability for bin i from the experimental data and <7j is the 
corresponding error. P- hcovy (z, a 2 , M 4 , M B ) are computed from Eq. 2.1, 2.12, 2.13. We assume 
the IH, so «2 is the only relevant unknown active neutrino parameter, and set the PMNS 
matrix to be Eq. 2.9. We include 8 bins each for LSND, MiniBooNE neutrino and MiniBooNE 
anti-neutrino data. We do not use MiniBooNE data for E < 475 MeV (L/E > 1.13 m/MeV) 
because of the unexplained excess in the low energy bins 9 . For the MiniBooNE anti- neutrino 
search, we performed two separate fits, using just the 2010 results or including the 2011 



8 The discrete sign in R can be absorbed by redefining the sign of a mass eigenstate field. For similar reasons 
it suffices to restrict < 02 < 7r. 



!) 



Comparison between data with E > 475 MeV and E > 300 MeV are performed in [56, 57]. 
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data. Before showing the results, we want to emphasize the differences between our fit and 
the global fits in [10, 11]. Though our fit is crude in the sense of not including all of the 
oscillation data available, it is a more direct test of the minimal mini-seesaw mechanism as an 
explanation for the short baseline oscillation anomalies. The global fits employ all physical 
parameters: the sterile neutrino masses, mixing angles and CP-violating phase, while the 
parameters in the minimal mini-seesaw are more restricted. Our fits use 5 real parameters 
characterizing the mini-seesaw mechanism: the complex angle z, active neutrino CP- violating 
phase «2 ; and two sterile neutrino masses. The mixing angles between the active and sterile 
neutrinos and the sterile CP- violating phase can be determined from these 5 parameters using 
Eq. 2.12. Even though we do not employ the full data set, we require the mixings between 
sterile and active neutrinos to be smaller than 0.15 in the numerical evaluation of the best fit 
to avoid conflicting with the reactor data or with other neutrino oscillation experiments with 
null results. 

The best fit results are presented in Table 2. They are in the right ballpark compared 
to the global fit values. In particular, the predicted parameters are quite close to the GL 
fit values. We also show the allowed regions of z = re lS at 68% and 95% C.L. from our 
simple fits in Figure 1, fixing the other three parameters to the best fit values. (The best fit 
constraint \U a i\ < 0.15 is not enforced in these contours.) The two fits using 2010 and 2011 
data give similar results, though the x 2 is larger for the 2011 data. The predicted transition 
probabilities as a function of L/E for both LSND and MiniBooNE at our best fit point are 
presented in Figure 2. Further details are given in the Appendix. 
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0.13 
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Table 2: Best fit points using the minimal mini-seesaw (MMS) formalism derived in the text. The 
first row uses the 2010 MiniBooNE anti-neutrino results while the second one includes the 2011 data. 
For comparison, we also show the global fit results KMS [10] and GL [11]. 



We have emphasized that one typically expects Mt to be non-zero and, in the mini-seesaw 
limit, to yield contributions to the active neutrino masses comparable to those induced by 
the mixing. Mt ^ of course reduces the predictabilty of the theory. Assuming that all 
of the mass terms are of the order of magnitude suggested by (2.5) this additional freedom 
could allow a NH or degenerate scenario for the active neutrinos. For more general scenarios, 
e.g., with different S fields or different dimensions of the effective operators for different 
families, one could allow one or more sterile neutrinos considerably heavier than an eV. The 
mixing of the light active neutrinos with the heavier sterile states would typically be strongly 
suppressed, but relatively large mixing angles (so as to still be relevant to LSND/MiniBooNE) 
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Figure 1: The allowed regions of z — re %e at 68% (red, dotted) and 95% (blue, solid) C.L., fixing 
ct2, M4, M 5 to the best fit values. Left: fit using the 20f0 MiniBooNE anti-neutrino data. Right: fit 
including 2011 data. They are almost identical. The best fit points are denoted by diamonds. 
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Figure 2: Predicted transition probabilities for LSND, MiniBooNE neutrino and MiniBooNE anti- 
neutrino oscillations from the KMS [10] and GL [11] global best fit points and from the best fit point 
using the minimal mini-seesaw formula. 
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would be possible by fine-tuning. 

2.5 Constraints and Implications 
2.5.1 Laboratory Constraints 

In this section we comment on two laboratory constraints on neutrino masses from tritium /3 
decay, 3 H — )- 3 Hee~z/ e , and neutrinoless double (3 decay {P(5q v ), i.e., nn — > ppe~e~. The shape 
of the electron spectrum near the endpoint in tritium j3 decay puts a stringent kinematic limit 
on the effective v e mass-squared 

m 2 p = J2\Aei\ 2 ml (2.15) 

i 

which is bounded to be in the range m 2 ^ = (1.1 ± 2.4) eV 2 from experiments in Mainz [58] 
and Troitsk [59]. The Karlsruhe Katrin experiment should improve the sensitivity down to 
around (0.2 eV) 2 at the 90% level [60]. The sterile neutrino state could have two effects on 
the electron spectrum. It could produce a kink in the spectrum of order \U e i\ 2 at energy 
EQ — rrii, where Eq = 18.6 keV is the endpoint energy and rrii is the sterile neutrino mass. As 
a result, it would be accompanied by a suppression of events at the endpoint of size 1 — |C/ e j| 2 . 
A detailed study has been performed in [61]. While current experiments do not exclude an 
eV-scale sterile neutrino with mixing angle ~ 0.1, Katrin would start to be sensitive to this 
parameter region. As shown in Table 3 the value of mp for our best fit point in the minimal 
mini-seesaw is 0.18 eV. 





3+2 (eV) 


3+2 MSS (eV) 


EXP (eV) 


S 


~ 0.2 
0-0.08 

~ 2 


0.18 


2.4 


(1 -2) -> 0.2 
(0.2-0.7) -> (0.01-0.03) 
(0.5- 1) -> (0.05-0.1) 



Table 3: Major constraints on the sterile neutrino masses. The second column lists typical values 
for generic (non-degenerate) 3+2 schemes, while the third are the results at the best fit point of the 
minimal mini-seesaw, mp and E are dominated by the effects of the eV-scale sterile neutrinos. The 
last column is the range probed by current and future experiments. 

The process f3/3o u constrains the effective Majorana mass in the presence of mixing be- 
tween light Majorana neutrinos, 

mpp = ^2(Aei) 2 mi = (M T )n, (2.16) 

i 

where the second equality is only true when all of the mass eigenstates are light compared 
to an MeV. In the particular minimal case we discussed above, Mt = and one expects 
m /3/3 = due to a cancellation between the light and eV-scale states. This is in contrast to 
generic 3 + 2 models [43], including the mini-seesaw with Mt ^ 0, in which mpp can vary 
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from to around 0.08 eV or 0.04 eV for the IH or the NH, respectively (both cases assuming a 
massless lightest state). The best current upper limit 10 is from the Cuoricino experiment [63] 
in the Gran Sassso Laboratory, which obtains mpp < (0.19 — 0.68) eV at 2a, with the range 
due to the nuclear matrix element. Future experiments should be sensitive down to (0.01 - 
0.03) eV [64, 65, 66, 67, 68, 69], which start to constrain M T . 

Other existing and future implications of mixing with eV-scale sterile neutrinos are con- 
sidered in [61, 70, 71, 72, 73, 74, 75]. 

2.5.2 Cosmological Constraints 

Though the 3 active + 2 eV-scale sterile neutrino model fits nicely into the mini-seesaw 
paradigm and could be explained easily by an effective theory containing a SM singlet with a 
TeV-scale VEV, it is nontrivial to accomodate two eV-scale sterile neutrinos into a consistent 
cosmological history. On one hand, current cosmic microwave background (CMB) and large- 
scale structure (LSS) observations show a slight preference for additional relativistic degrees 
of freedom beyond the SM prediction. The combination of the CMB, measurement of the 
Hubble parameter from HST observations and baryon acoustic oscillation (BAO) data shows 
that the effective number of relativistic degrees of freedom N e f f is greater than the SM value 
3.046 at 98.4% C.L. [76]. By further combining with the Atacama Cosmology Telescope 
(ACT) power spectrum measurement, the number is estimated to be N e ff = 5.3 ± 1.3 (68% 
C.L.) [77]. Similar analyses that show preferences for light sterile neutrinos can be found 
in [78, 79]. However, it is argued in [80] that once the prior-dependence in the cosmological 
analysis is removed, the latest cosmological data show no evidence for deviations from the 
SM value of neutrino species. In the near future, data from the Planck satellite will be able 
to reduce the error in N e ff to 0.20 and settle the issue [81]. 

On the other hand, the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data, 
small-scale CMB observations, Sloan Digital Sky Survey (SDSS) 7th data and the present 
Hubble constant, set a 95% C.L. limit on the sum of two light sterile neutrino masses [13] 

5 

£ 5 = ^2\rm\ < 0.9eV (95%C.L.), (2.17) 

which implies S = ^ i |mi| < 1 eV for the non-degenerate IH, which is small compared to the 
typical values in Table 3. 

It is also well known that if two eV-scale steriles are fully thermalized prior to neutrino 
decoupling, they would increase the Hubble rate and result in an overproduction of 4 He in the 
epoch of big bang nucleosynthesis (BBN). A combined analysis of 4 He and deuterium data 
sets a 95% C.L. limit on N s < 1.26 [82] with a best-fit value N s = 0.86. So while one fully 
thermalized neutrino is slightly favored over none, two fully thermalized sterile neutrinos with 
masses around 1 eV are not fully compatible with all the astrophysical data unless certain 
assumptions of these analysis are relaxed. 

10 One experiment [62] claims to observe a positive signal, corresponding to mpp ~ 0.3 eV, but this has not 
been confirmed. 
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There are several (highly speculative/creative) loopholes to these cosmological constraints. 
One possibility is that there exists a common non-zero chemical potential for the active neutri- 
nos £ [83, 84, 85]. In the BBN epoch, a non-zero chemical potential would drive the reactions 
n + v e <->■ p + e~ , n + e + <->■ p + z7 e in one direction and change the ratio of neutrons to protons 
in thermal equilibrium to 



The additional suppression would compensate the effects of an increase of the relativistic de- 
grees of freedom. Assuming two light sterile neutrinos, a positive non-zero chemical potential 
is needed 0.03 < f < 0.14 at 95% C.L. with a best-fit value of £ = 0.064 [82]. Such a large 
chemical potential would correspond to an enormous lepton asymmetry of the same order 
as £ (~ 10 s x the baryon asymmetry), but might be generated from a much smaller initial 
asymmetry due to non- linear effects associated with the sterile neutrinos [86, 87, 88]. 

Another possibility discussed in the literature is that the sterile neutrino has a time- 
dependent mass [89, 90, 91]. This could be achieved if the singlet that generates the sterile 
neutrino mass has a medium-dependent VEV, i.e., a VEV proportional to the number density 
n v of the active neutrinos: (S) oc n u . Then M% = Mg{l + z) n in our effective model, where 
Mg is the mass in the epoch with redshift z and n = 6 if the sterile Majorana mass is 
generated from a non-renormalizable operator S 2 N 2 11 . Then at the recombination epoch 
with z ~ 1100, the sterile neutrino mass becomes 10 9 GeV for the present mass of 1 eV. 
Thus the production of sterile neutrinos are highly suppressed in the early Universe and their 
abundance is negligible. 

A third possibility is to have late time phase transitions to suppress the sterile neutrino 
masses and mixings until after they decouple at T < 1 MeV. This could be achieved in models 
with low-scale breaking of global symmetries [92]. 

3. Heavy Sterile Neutrino 

3.1 The Allowed Parameter Region 

There is another interesting possibility [14], that the LSND and MiniBooNE results could 
be consistent with the paradigm with one eV-scale sterile neutrino and another very heavy 
sterile with mass ranging from O(10) eV to O(10) GeV. In the presence of a coherent heavy 
sterile neutrino 1^5, x§ varies very rapidly and must be averaged over. Furthermore, the 
heavy neutrino could be produced incoherently and/or with a reduced phase space. If it is 

11 The discussion here is slightly different from that in [45], in which n = 3, as the sterile mass is assumed to 
be from a renormalizable operator SN 2 . It is also assumed in [45] that the Dirac mass mixing the sterile and 
active neutrinos is generated by the Higgs VEV and does not change with time. Thus the induced contribution 
to the light mass MoMg 1 and the active-sterile mixing angle MoMg 1 decrease at early times. However in 
our effective theory, the Dirac mass also decreases with time Mn oc (1 + z) 3 and thus the induced contribution 
to the light neutrino mass is constant while the mixing angle still decreases with time. 




(2.18) 



- 13 - 



sufficiently heavy, above the energy of the experiment, it is not produced at all. Following [14], 
we introduce two parameters: the mixing ratio r and the CP odd quantity /?, defined as 

\u e5 u; 5 + u e ,u;,\ 



r = 



\U e4 U, 



i' i 



so that 



R = 1 tan" 1 ( smdpesU^l 

2 Vl^e4^4|+COS(5|C/ e 5^5| 



D 2i/3 _ U esU* 5 + LV e4 LT* 4 



(3.1) 



re^ = g ~ (3.2) 

Also taking into account the possibly reduced phase space by including another free param- 
eter a, which varies from (no production) to 1 (no suppression), we have the appearance 
probability 

*V* e = \Uei\ 2 \U^\ 2 ((1 - r) 2 + a ((1 - r) 2 + 4r sin 2 $) + 4r sin 2 (x 4 + /3)) , (3.3) 

with /3 — > —f3 for anti-neutrinos. 

A x 2 fit h as already been performed by Nelson [14], which indicates that the modified 3+1 
scheme with CP violation could possibly account for the short baseline anomalies. Here, we 
perform a similar fit to allow easy comparison with the minimal mini-seesaw fits in Section 2 
and to incorporate the MiniBooNE 2011 data. We follow a similar procedure to [14], with 
the probability given by Eq. 3.3. We numerically evaluated a \ 2 function similar to Eq. 2.14 
with the free parameters \U e 4U^\,r, (3,m 4 i while fixing a. (The fit does not vary much with 
different values of a.) We include 8 bins each from LSND and MiniBooNE neutrino and 
anti-neutrino data as in Section 2.4. Following [14], we also include one bin for KARMEN 
and one for NOMAD [93], with the experimental errors chosen to be the 90% upper limit on 
the average oscillation probability. We also impose the constraint \U e4 U^\ < 0.03 to avoid 
too large effects in reactor and other experiments. The total number of fit points included is 
26 with 4 free parameters. For the MiniBooNE 2010 data, at a = 0, the best fit point has 
PeiU^l = 0.03, r = 1.10, = -0.15, Am^ = 0.49 eV 2 with a x 2 = 21.5. For a = 1, one finds 
\U e4 U^\ = 0.03, r = 1.04, /3 = -0.14, Am^ = 0.48eV 2 with a x 2 = 21.9. The results, which 
are consistent with those in [14], are presented in Figure 3, which show the x 2 an d | £44^4 1 
distributions as a function of Am| 1; with \U e iU^i\,r, (3 chosen to minimize x 2 a t each Am^. 
We also show \U e bU^\, as calculated from Eq. 3.2. We caution that, similar to Section 2, this 
fit does not incorporate all data, except for the imposed upper limits on the mixings, and 
should only be viewed as evidence for the modified 3+1 scheme with CP violation as a solution 
to short baseline neutrino anomalies. The probability distributions for the short baseline 
experiments at one particular point in the parameter space with a small x 2 are shown in 
Figure 4. The quality of the fit gets worse when including 2011 MiniBooNE data. For instance, 
at a = 0, the best fit point has \U ei U^\ = 0.03, r = 1.08, /3 = -0.12, Am 2 ^ = 0.47 eV 2 with a 
X 2 = 27.5. 
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Figure 3: The minimal values of % 2 and best fit values of IL^E/^I and \U e ^U^\ as a function of 
Am^ in the scenario with a light and a heavy sterile neutrino. The constant value of | C/ e 4 1 for 
small A 771^ is due to the imposed constraint IL^C/^I < 0.03. 
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Figure 4: The predicted probabilities for LSND, MiniBooNE neutrino, and MiniBooNE anti-neutrino 
oscillations in the Modified 3 + 1 scheme at Am^ = 0.6eV 2 , a = 0, and | t/ e 4 1 ' r > P chosen to 
minimize \ 2 - Also shown for comparison are the global best fit results in the 3 + 2 model. 



3.2 Constraints 

One crucial feature of the results in the previous section is that the heavy sterile neutrino, 
independent of its mass, needs to have large mixing angles with v e and z/^, e.g., \V e ^U^\ ^ 
1CP 2 , to yield large enough CP violation required for the explanations of the short baseline 
anomalies. Similar to light sterile neutrinos, a heavy sterile neutrino with large mixing is 
also constrained by CMB, LSS and BBN data. There are also several other laboratory and 
cosmological constraints. All of the constraints are summarized in Figures 1 and 2 in [45, 94], 
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to which we refer the reader for more details. Below we just review the most stringent ones 
in different sterile neutrino mass ranges. 

For a sterile neutrino with mass above 1 eV but below 5 keV, LSS sets the strongest bound 
on the energy density stored in the sterile neutrino as a function of its mass, and rules out 
mixing of order 0.1. The heavy sterile neutrino could decay to an active neutrino and a photon 
at the one-loop order: v s -4 v a + 7 with ~ ms/2. Thus the X-ray emission from galaxy 
clusters such as the Virgo cluster or from all possible sources seen as the diffuse extragalactic 
background impose a strong limit on the flux from the decay. The bound is strongest in the 
mass range from a few keV to around a few hundred keV, restricting the mixing-squared to be 
smaller than at least 10~ 5 . If the light active neutrinos from the decay are produced between 
the active neutrino and photon decoupling time, they would contribute to the energy density 
of the relativistic species and modify the CMB angular power spectrum. The CMB spectrum 
rules out a heavy sterile neutrino with mass from 0(100 keV - 10 MeV) with mixing > 10~ 2 . 
Sterile neutrinos with mass 1 MeV < m, < 400 MeV are constrained by the non-observation 
of additional lines in the charged lepton spectrum and the decay products from a heavy 
sterile neutrino produced by meson decays such as (tt + , K + ) — > fi + u^ [95], which also rule out 
mixing of order 0.1. For an even heavier sterile neutrino with mass > 0.1 GeV, there are weak 
constraints on the mixing < (0.03 — 0.1) depending on the mass from flavor changing neutral 
currents and violation of lepton universality [96, 97]. In addition, LEP searches for a SM 
singlet neutral heavy lepton Ni in the channel Z — > Nivi exclude the mixing-squared down to 
below 10~ 2 in the mass range between 1 GeV and 90 GeV [98]. Thus, without modifying any 
cosmological assumption, there is only room for a heavy sterile neutrino with mixing ~ 0.1 
for mass in the range (0.4 - 1) GeV . 

Finally we comment that there are two additional ways to avoid the cosmic constraints 
for heavy sterile neutrinos. The first scenario has a low reheating temperature after inflation, 
Tr <C 100 MeV to dilute the sterile neutrino density and relax the cosmic bounds, but Tr > 5 
MeV to preserve a successful BBN [99]. Another possibility is to have the sterile neutrino 
decay to a massless or very light (pseudo-)Goldstone boson, e.g., Majoron, with a very short 
lifetime t«1s [100]. 

4. Conclusion 

Motivated by short baseline neutrino oscillation anomalies, we study two scenarios with light 
eV-scale sterile neutrino(s). The first has two eV-scale sterile neutrinos and a mass pattern 
arising naturally from the mini-seesaw mechanism. The mini-seesaw could be generated 
by a set of higher-dimensional operators involving a new TeV-scale VEV, which could be 
related to TeV-scale flavor models. In the limit where the active neutrino Majorana mass 
vanishes, a simple formula is derived relating the 6 complex mixing parameters between the 
sterile and active neutrinos to only 5 real parameters. A simple numerical fit is performed 
based on this minimal mini-seesaw formula, indicating that it could account for the short 
baseline anomalies. We also review the current laboratory and cosmological constraints and 
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implications for this scenario. The second paradigm has one eV-scale sterile neutrino with 
a CP-violating phase generated by integrating out a heavy sterile neutrino. However, for 
large enough CP violation to explain the anomalies, the heavy sterile neutrino has to have a 
large (fine-tuned) mixing with the active ones, which is in strong tension with astrophysical 
and laboratory observations. Ongoing and future neutrino oscillation experiments, as well 
as cosmological observations such as the Planck mission and laboratory searches for neutrino 
mass effects in [3 decay and neutrinoless double /3 decay will hopefully confirm or exclude the 
existence of light sterile neutrinos as an explanation of the short baseline anomalies. 
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A. More details of the Fit with the Minimal Mini-Seesaw Parameterization 

In Section 2.4, we present a simple fit based on Eq. 2.12, using Mathematica to calculate 
the x 2 an d find the best fit point. The result, listed in Table 2, uses data from LSND and 
MiniBooNE, as well as a constraint that the active-sterile mixing matrix elements should not 
exceed 0.15 in magnitude. If we had not required \U a i\ < 0.15, some of the mixings would 
have been even larger at the best fit point (the largest was \U e ^\ ~ 0.3). This is a reflection 
of the well-known fact that there is significant tension between the short baseline data and 
other reactor and accelerator constraints [9, 10, 11, 56, 57]. Nevertheless, the restricted fit 
has a quite reasonable x 2 > as seen in Table 2. 

We compare the contours of mixings and of the allowed region (not imposing the mixing 
contraint) in the (r, 9) plane at 95% C.L. in Figure 5. It shows that most of the allowed region 
leads to reasonably small mixings < 0.2. We also show the dependence of x 2 as a function of 
the phase a 2 in the active sector in Figure 6. The x 2 does not change much over the entire 
range, with two minima at 012 = 1.2, 2, where the latter has a slightly smaller x 2 - 12 
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